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the framework of the model, factorization in pseudoscalar transitions holds 
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I. INTRODUCTION 



The factorization approximation is widely used in the heavy quark physics. As a calcu- 
lational scheme, it is very convenient and there are numerous models based on it. However, 
the theoretical basis for the approximation itself is not completely understood. The usual 
motivation for it is a color transparency argument due to Bjorken |Q which suggests that 
the quark pair produced in the decay does not have enough time to evolve to the real-size 
hadronic entity, being a small size bound state with a small chromomagnetic moment. That 
is taken to imply that the QCD interaction is suppressed. For the heavy-to-heavy-light 
transition (we shall hereafter refer to the color allowed decay B — > D^n^ for the sake of 
clarity) the matrix elements in question are 

M{B{pb) ^ D{pD)7c{q)) = {TcDlnmB) = ^V,kV,^{[C, + ^C2)Mi + 2C2Ms} = 
%VcbVdu{{Ci + ^C2)(7rZ}|c(0)7^(l +75)&(0) (1(0)7^(1 + 75)^(0) + 

2C2{nD\c{0)^^{l + ^,)t%{0) diO)ril+l5)taUmB)} (1) 

Here the usual notations are used 

n = CiO, + C202, {nD\ O, \B) = Ml, {nD\02\B) = ^^Mi+2M8 
Oi = Q7^(l + 75)6^ dkYi^ + 75)w' and O2 = Q7^(l + 75)&' 47^(1 + 75)^' (2) 

where we have explicitly written out color indices. In the usual factorization approximation 
one simply separates out currents by inserting the vacuum state and disregards any QCD 
interactions between them on the basis of Bjorken's conjecture. That implies that the 
A4s contribution in (|l]) is dropped, as it has no factorizable color singlet form. New non- 
factorizable corrections appearing in the A^g matrix element can arise if a hard or soft gluon 
is emitted or absorbed by the escaping pair. Obviously, a more quantitative estimate of 
these contributions is needed. 

The first attempt to put a factorization approximation on the solid theoretical base was 
made in . The authors considered a limit of large energy transferred to light quarks 
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as rrib — >■ oo. Here v is the velocity of the decaying heavy quark pb = rubV + k and q is the 
momentum transferred to the hght quark pair. In this hmit the energy scales as a heavy 
quark mass. Proposing a Large Energy Effective Theory (LEET) as an extension of the 
Heavy Quark Effective Theory and making use of a convenient gauge condition n ■ A = 
(with n being a null vector, n'^=0) they proved factorization for the physical amplitude 
dominated by collinear quarks. It is not clear, however, how the pair of coUinear quarks 
hadronizes into the pion. Also, the validity of the LEET itself is not yet established. 

The important step towards the sought theoretical description was made in 0. The au- 
thors considered a specific small velocity limit where m^ — ~ const which implies different 
scaling of (Q): ~ const as — * oo. It turns out that this limit is theoretically "clean" 
for the application of the QCD sum rule method and provides a theoretical justification 
of the "rule of discarding l/N^ corrections" 0] on dynamical grounds. Unfortunately, the 
actual B decays are at the borderline of the formalism, bringing large uncertainties into the 
estimation of the matrix elements. 

In what follows we develop a QCD-based model which works in the limit E ~ rub as 
nib oo. The idea is to study the propagation of the light quarks emerging to form the 
final hadron while passing through the (perturbative or non-perturbative) gluonic fields 
present in the decay. Interestingly, results can be obtained which are independent of the 
nature of the gluonic fields. The light quark pair, produced at x = 0, interacts strongly 
with the background gluonic field while escaping from the point they were produced and 
hadronizes into the pion. Since the energy transferred to the quark pair is large it is not 
necessary that the hadronization occurs at a; ~ 0. As a result, it is not sufficient to take 
the operators of the lowest dimension in the Operator Product Expansion (OPE) but rather 
one must sum up a whole series of operators. The novel technical point of our method is 
to use introduce "generating (distribution) functions" which incorporate this infinite series 
of matrix elements, and then to use Heavy Quark Symmetry to restrict the form of the 



matrix elements enough that useful statements can be made. In principle, these functions 
could be modeled or perhaps fixed by other measurements. However, at the leading order in 
= 0, ~ the color octet matrix element vanishes due to cancellations and factorization 
holds. The method should be able to be extended beyond leading order, allowing one to 
also see corrections to factorization. 

II. THE OCTET CONTRIBUTION 

The amplitude for the singlet or octet part is 

Mi,s = {D7r\c{0)^^{l + 75)T6(0) Jt'^II + l5)Tu\B) 

= {D\ciOh,{l + ^,)Tb{0)V,\B) (4) 

with T = being the appropriate color matrix. The quantity describes the tran- 

sition of the light quarks into a pion in the presence of the gluonic fields which are either 
internal to the B,D hadrons or produced during the transition. The idea is to treat these 
gluonic fields as external fields with respect to the light quark pair produced in the decay 
of virtual H^- meson. From the diagrams Fig.l one may write down the expression for one 
gluon matrix element of the produced quark pair in this general background field. This gives 
a contribution to the octet matrix element only: 

= /rf^x(7r|J(x)A"(x)4^^(x, 0)7^.(1 + 75)t^.M(0)|0)- 

{tgs) I f/^x(7r|rf(0)7^(l + l^K.SpiO, x)JA\x)tiu{x)\0) (5) 

where we set = A^ = the external field. As one can see, in this approximation we neglect 
higher Fock state contributions. Our normalization is such that Tr f^t^ = 6°'^/2. Thus we 
can rewrite the equation as 

V, = -'fJ d^x(Tr^-{x)SF{x, 0)7^(1 + 75)G(0, x) + 

Tr7^(l + 75)^f(0, x)^'^{x)G{x, 0)) (6) 
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with non-local correlators 

G(0,x) = {■K\Tu{Q)d{x)\Q) 

G'(x,0) = (7r|Tn(a;)J(0)|0) (7) 

Shifting the x-dependence of the second correlator 

(7r|Tn(a;)rf(0)|0) = (7r|Te'^"u(0)e-*^X0)e*^"e-*^"|0) = e''?"G(0, -x) . (8) 

We shall work in the Fock-Shwinger gauge which gives us the advantage of expressing 
the Afj^ in terms of gauge-invariant quantities (such as the gluon stress tensor G^^): 

x^-A^' = 0, = daa XpG^^(ax) . (9) 

We treat the light quarks as massless, so the bare propagator of a light quark is given by 

^'^^-y^ = irH{x^-yw 

This allows us to rewrite as 

-zV, = -iVi ^-iV2,= ^-^!d^x ada [^G%{ax)Trastl^.{l + 75)^(0, x)- 

2^G^„(ax)Tr7^(l + 75)cT.nG(0, -x)e*^-'] (11) 

The correlator (7(0, ±x) entering the expression above is non-perturbative and nonlocal but 
can be related to the wave function of the forming pion. Following [^] we expand it near the 
light cone (x^ ~ 0). The general form of expansion is 

Tu{Q)d{x) = Y.Cn{x^ - ie)V^x^K..x^-{dV''D^,...D^^u) + 

n 

Y.CUx^-ie)T^x^\..x^-{dT^D^,...D^^A^^^,...A^^u) (12) 



n,m 



with Fq being the full set of Dirac matrices. We would like to note, however, that since 
(7(0, x) (or G{x, 0)) contributes to the trace, only 7^ and 757^ terms enter the final expression 
for (7(0, x) {G{x, 0)). In Eq. ([T^), the non-gauge invariant terms are exactly zero by the FS 
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gauge condition (Bl) which allows us to use analysis similar to B. Moreover, in this gauge 



xi'D^^ = x^d^^ in (|T2|). Thus, 



n 

Parameterizing unknown matrix elements in the conventional form 



(13) 



(7r|rf ™^,...D^„m|0) = i?„g%,...g^„ (14) 
where g is a momentum of the pion we arrive at the simple form for the G-correlators 



(15) 



(7r|rM(0)rf(x)|0) = (r„g°) j dze"^''(j){z) 
{7r\Tu{0)di-x)\0) = (r„g") / dze'^'-''^'"'(j){z) 



with the function (j){z) introduced through its moments 



dzz'^Mz) = (-iY'n\C„B„ 



(16) 



The zeroth moment gives the overall normalization of (f){z): J dz (f){z) = /^r/ {2^/Nc) with f^^ 
being the decay constant of vr-meson. Inserting (|15]) into the (|^) we obtain 



Qs f , ,A 0{Z) 



^V,, = 5 / dzd^xr^^-^ 



2fi — 2 



dzd X- 



x^ 



^'-'^'"^{x.q ■ A\x) - 4(x)(x ■ q) + zest,aAl{x)xtqa} (17) 



This gives for the octet amplitude of i? — Dti transition: 

M'i^ = Ml + M2 = (D|c(0)7m(1 +75)ta[(vr|J7^(l +75)w|0)^]6(0)|5) = 
x^^q,. - (x ■ q)g^^ - ie^^^pXaqp {D\c{Q)T ^A^{x)h{Q)\B) - 
x^qu - {x ■ q)g^,u + ie^uapXaqp {D\c{Q)T f,Ay{x)h{Q)\B) 



i9s 



dzd^x^e'''"' 



x^ 



dzd X- 



TC^ J X^ 

where f = 7^(1 + 75) and fA^ix) = A{x)^. 



:i8) 



^Note that the transferred momentum q can be written as g = En where n is a null vector (as 
in the approximation of collinear u and d - quarks) and E is energy transferred to the light quark 
system (it is large in the kinematical limit under consideration). Observe that contrary to Q there 
are terms in (^) which are not proportional to n • A. 



This amplitude is also nonlocal, so we use the same idea as in Eq. (jT^) expanding it 
around ~ and introducing "generating functions". It is clear that the only matrix 
element to be parameterized is 

{D\c{0)T,A,{x)b{0)\B) (19) 
In FS gauge one can expand A^{x) about x = using "gauge invariant" decomposition 

OO -| 

M^) = E (;^:^^"-^^"^''K-^^.G,,(o)) (20) 

This implies the matrix elements 

{D\c{0)T^A,{x)b{0)\B) = 
2,. ^n=o j^x-^^...x-^-xP{DmT^D,^...D,^\D^,D,]bmB) = 



2g, 



J2n=0 (n+2)n!"^ ^---X X'^ M ^y-^^^^^y^p^ (21) 



In order to parameterize the matrix elements appearing in the expression above we note 
that there are a total of n + 3 indices and that the whole expression is anti-symmetric with 
respect to p ^ u interchange. Also, there are only two vectors (or their combinations) 
available, so we choose them to be pb and q. This leads to 

EILo biqpPBu,--PBu,(lu,+,--qu„ {pBpqu - PBuQp) (22) 

with Oj and hi being unknown coefficients. Performing the contraction and dropping terms 
proportional to q^ and we have for the first matrix element in ( |l8l) 



Ml = -2^/^'a;rfz^e--''E;r=oEr=o (;i^{-2a.(pBg)(PBx)(xg)- 

bi{xqf{pBq) + aim\{xqY + hi{xqY{pBq)]{pBx)\xqY-''- (23) 

Note that in this limit the dependence of the matrix element on hi cancels out. These sums 
can be parameterized in terms of some distribution function f{n,v). Using the binomial 
formula, we obtain 



(24) 



where C\ is a binomial coefficient. This gives for the matrix element 

Ml -- 



27r 



Performing resummation with respect to n we arrive at 



n + 2)n! 



/ i/) 



(-0" 

n+2 



/(^, '^) 



(26) 



This results in 



Ml 



27r2 



dzduda(j){z)'ip{a, v) m%q^,q^ - 2{pBq)PB^qu 



,4 X^Xt, Maupg+{a+z)q)a 



d X e 



(27) 



Similarly, for the second matrix element one obtains 
1 



Mo 



27r2 



d^dz/cia0(z)^(a, I/) [mlq^q, - 2{pBq)pB^,qu\ J (28) 



At this stage we have classified the infinite number of the unknown matrix elements in terms 
of a few continuous functions. The strategy at this point in general would be to consider 
the structure of these functions and make reasonable smooth guesses as to their form. For 
example (j){z) is related to wavefunction used in the studies of the perturbative limit of the 
pion form factor and has the asymptotic form 0(z) = 2/^2;(l — z)/^/6 0. However, we find 
that an important cancellation occurs which is independent of the form of these structure 
functions. Taking a/3 = (a + z)qi3 + aupBjj, &/3 = (1 + + a)g/3 + aupBp and integrating over 
X we obtain for the Als^'':^ 



^The key integral is 



a X- 
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Ai^s^ = —i J dvdail)i{a, v) J dz(f){z) 
{j^[4:{pBq){PBa){qa) - 2m%{qaf - 2a^{pBqf 
w[^{pBq){pBh){qh) - 2mUqbf - 2b\pBq)^]} = (29) 

at the leading order in for any and ip. This imphes that the semi-soft-one-gluon matrix 
elements do not bring new contributions to the factorization approximation result in the 
kinematical limit under consideration. 

We would like to point out that the two-pseudoscalar final state is somewhat special: for 
instance, B —>■ Dp does not fall into the same category - there is a good chance that under 
the similar conditions there is no complete cancellation of the diagrams. This is due to the 
fact that there can be extra terms in ( ^2]) that include polarization vectors of p meson. 

In conclusion, we have constructed a QCD-based model for the factorization approxima- 
tion in the limit E ~ ^bmt^oo- the framework of the model the complete cancellation 
of the non-factorizable terms is observed at the leading order in g^, x^. The OPE is re- 
summed by introducing two generating (distribution) functions that can be independently 
parameterized. The model gives a consistent way to estimate non-leading (non-factorizable) 
contributions. 
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